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Abstract. We study a one-dimensional particles system, in the overdamped limit, 
where nearest particles attract with a force inversely proportional to a power a of their 
distance and coalesce upon encounter. The detailed shape of the distribution function 
for the gap between neighbouring particles serves to discriminate between different laws 
of attraction. We develop an exact Fokker-Planck approach for the infinite hierarchy of 
distribution functions for multiple adjacent gaps and solve it exactly, at the mean-field 
level, where correlations are ignored. The crucial role of correlations and their effect 
on the gap distribution function is explored both numerically and analytically. Finally, 
we analyse a random input of particles, which results in a stationary state where the 
effect of correlations is largely diminished. 
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1. Introduction 

In a recent paper [1], we have mapped the late-time evolution of the conserved 
Kuramoto-Sivashinsky (CKS) equation [2], 



used to describe, e.g., the dynamics of meandering steps in crystal growth or that of 
wind driven sand dunes, to a reacting particles system. The long-time profile u{x,t) is 
made up of pieces of a universal parabola (solution of the equation u + Uxx + = c) , 
joining in regions of vanishing size and diverging curvature whose locations define the 
profile without ambiguity, and that one may associate with 'particles'. The projected 
particles on the x-axis form an overdamped one-dimensional system; nearest particles 
attract one another with a force inversely proportional to their distance, and particles 
coalesce irreversibly upon encounter. This process gives rise to a reduction in the particle 
density, i.e., to a coarsening process. 

In this paper, we generalize the above reaction model to the case where the force 
between nearest particles, separated by a distance i, is proportional to (a = 1 
corresponding to the CKS equation). Reacting particles systems such as this, evolving 
deterministically, occur often enough but have been studied sparingly. Coarsening in 
spin systems has been studied for exponential forces [3] between the kinks and for power- 
law forces (between all particle pairs) [4, 5]. The deterministic Kardar-Parisi-Zhang 
equation corresponds to an interacting particles system [6], and the (non conserved) KS 
equation was shown to be equivalent to coalescing particles traveling ballistically [7], as 
is also the case for the Burgers equation at high Reynolds numbers (shocks representing 
particles) [8]. Similar ballistic systems have been studied by Redner et al. [9]. Simple 
deterministic reaction schemes, such as the "cut-in-half model, have been devised by 
Derrida et al. [10], but exact results are few even in those ideal cases. 

The aim of this paper is to give a thorough description of these models with 
varying a. The time-dependence of coarsening can be easily found using dimensional 
analysis, however, it provides too weak a characterization of the kinetics (e.g., diffusion- 
limited one-species coalescence coarsens in the same way as our model for a — 1). 
We therefore focus on the distribution function g{i,t) for the gap i between nearest 
particles. We derive an exact hierarchy of Fokker-Planck equations for the probability 
density functions (pdf) {pm}m=iJ finding m consecutive gaps of given sizes at time t. 
Truncating this hierarchy at the mean- field level, ignoring correlations between adjacent 
gaps, we solve for the (approximate) single gap distribution function. The shape of this 
pdf is a much better discriminant between various models, and we focus in particular on 
the exponent P characterizing the small gap behavior, g{i) ~ The role of correlations 
is explored numerically and analytically, using the Fokker-Planck equations to derive 
exact constraints that the gap distributions must satisfy when correlations are fully 
taken into account. Finally, we test the mean-field approximation for the case where 
a homogeneous random input of particles produces a steady state that is considerably 
less affected by correlations. 




(1) 
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Apart from the intrinsic interest in deterministic reaction models, our particular 
system might also be relevant to the dynamics of vicinal surfaces undergoing step- 
bunching instability [11] — a process where parallel straight steps do not move uniformly 
but rather bunch together in groups separated by large terraces. In this picture, particles 
represent steps. One of the several ways that steps interact is through diffusion, giving 
rise to an effective attraction decaying as 1/(6 + £), where b is some characteristic 
length. For large terraces, the model a = 1 is recovered, while the opposite limit, £ <^ b, 
corresponds to our model with a — —1. 



2. The Model 

Consider an infinite system of particles on the line, located at {xi{t)}'^_^. The system is 
overdamped, and nearest particles attract one another with a force inversely proportional 
to a power a of the distance: 

dxi A ( \ 1 



- ■ (2) 

The prefactor A carries units of (length) "+^/ (time) (or L"+^/T, symbolically), and 1/a 
has been introduced to enable discussion of negative a, as well as the limit a — > 0. 
Prom (2) one finds that the gaps between particles, £i — Xj+i — Xi, obey 

d£^^Ar^_l j_\ 
dt a\eu J ' 

The particles coalesce upon encounter, or analogously a gap is removed from the system 
once it has shrunk to zero. Thus, the number of particles decreases and the typical gap 
between particles grows as the system coarsens with time. 



2.1. Linear Stability Analysis 

To analyse the system's stability let us perturb a uniform configuration, 

(^^ = {i)+e^, ir (^1 - (4) 

and substitute in (3) to yield 

A 



-I — + ~ 2ei) , (5) 

or, upon passing to the continuum limit (for large i and t), — > e{x, t), 

—e(x,t) — rY-;r-:^e{x,t) . (6) 

dt ^ ^ {iy+^dx^ ^ ^ ^ ^ 

Thus the perturbations obey a diffusion equation with negative coefficient, and the 
system is unstable regardless of the sign of a: the prefactor 1/a introduced in (2) 
ensures this outcome. For a > the attractive forces between particles separated by 
short gaps grow ever stronger, resulting in the eventual coalescence of the particles, 
and corresponding to singularity buildups in the diffusion equation. For a < 0, the 
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forces are repulsive and short gaps shrink into singularities due to the dominating forces 
of the larger surrounding gaps. In the demarcating limit of a ^ the forces depend 
logarithmically upon the distance: 

^ = Aln^i^^; ^ = Aln—^. (7) 
at Xi^i — Xi at 

2.2. Simulations 

Simulations have been performed by numerical integration of Eq. (2) for systems of 10'^ 
particles (for a > 0) and 10^ particles (a < 0). The calculations were carried out on a 
Linux AMD64 cluster comprising 32 processors. Data were obtained by averaging over 
thousands of independent runs. The long-time asymptotic limit is reached sooner for 
positive a, so typically the runs proceeded until the number of particles dwindled to 
500-700 for a > 0, or 50-90 for a < 0. 

3. Dimensional analysis 

Numerical integration of the equations of motion, starting from different initial 
conditions, suggests that the original state of the system has no effect on its long-time 
asymptotic behavior. In the long-time regime, therefore, the only dimensional physical 
parameters determining the typical gap £ (a length, L) are the elapsed time t (T) and 
the constant A (L"+^/T), so one expects 

£~ (8) 

This is indeed confirmed by numerical simulations. 

This scaling argument breaks down as a becomes smaller than —1. As a — > — 1 
coarsening proceeds ever faster and takes place at a finite time (even in an infinite 
system) in the limit a — —1, and the initial configuration becomes then important. 
Independence from initial conditions, however, is a prerequisite for the scaling behavior 
that wc assume (and find numerically) in order to analyse the gap distribution in 
Section 4. Accordingly, we limit the discussion to a > — 1. 

If in addition to the coalescence reaction particles are input randomly, at rate R 
per unit length per unit time (or 1/LT), the system achieves a steady state that is 
independent of its initial condition. The typical distance at the steady state is then 
dependent upon A and R alone, so 

^\ l/(2+a) 



Physically, the characteristic time for the shrinking of ig, il'^°'/A, is then balanced by 
the time l/isR for getting a new particle into the interval. For the important case of 
a — 1, the exponents i ~ i^/^ and 4 ~ R~^^^ agree exactly with the well known cases 
of diffusion-limited coalescence and diffusion-limited coalescence with input (where the 
particles diffuse freely without damping and without an attractive force) [12]. 
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4. Gap distributions and Fokker-Planck approach 

We see that the coarsening kinetics is not sufficient to distinguish between such different 
systems as diffusion-limited coalescence, annihilation, and our model with a = 1. A full 
characterization of a particles system on the line is provided by the hierarchy of joint 
pdfs for the multiple-gap correlation functions pm{h,^2, ■ ■ ■ ,£m,t) — the probability of 
finding m consecutive gaps of lengths £i, £2, • • • , -^m at time t. This is fully equivalent to 
the more traditional hierarchy of multiple-point correlation functions pk{xi,X2, ■ ■ ■ ,Xk,t) 
for finding k particles at Xi,X2, ■ ■ ■ ,Xk at time t [13]. In the following, we develop a 
Fokker-Planck approach for obtaining the pm- For practical purposes, however, the 
single gap distribution pi = p{i, t) provides already with enough detail to distinguish 
between various reaction models. For example, the tail of the distribution falls off 
as a gaussian, p ~ exp(— ci£^/t), for diffusion-limited coalescence, but as a stretched 
exponential, exp(— 02^^''^), for coalescence with input, while for small ^ one finds p 
with /3 = 1 in both these cases but /? ~ 1.3 for our model with a = 1. For other values 
of a one finds that P has a non trivial dependence upon a (Fig. 1). 
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Figure 1: The exponent (3, that describes the small gap behavior of the gap distribution 
function, p{i, t) ~ as a function of the force exponent a in our deterministic reaction 
model. The solid line indicates the mean-field prediction of Eq. (23), j3 = a. 
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Let nm{ii,i2, . . . , im, t)diidi2 ■ ■ ■ dim denote the number of m consecutive gaps of 
lengths between £.j and ij + dij (j = 1, 2, . . . , m), at time t. m = 1 corresponds to the 
single gap distribution nidi = n{i,t)di — the number of gaps between length £ and 
i + di, at time t. The total number of gaps (or particles), 

poo 

N{t) = / n{i, t)di^l/ir^ r^K^+o.) ^ (10) 
Jo 

decreases with time, since the total length of the system A ~ N{t)i is constant. Likewise, 
the normalization for m > 1 is 

/»oo poo 

/ •••/ nm{ii,...,im,t)dii---dim = N{t), (11) 

JQ Jo 

assuming periodic boundary conditions. For free boundary conditions the result is 
N{t) — m + 1, but (m — 1) can be neglected in the thermodynamic limit of N(0) — > oo. 
The joint probability density, p^, is related to n^n via the N{t) normalization: 

Pm(il, ...,im,t)^ N(t)-^nm{il, ...,im,t). (12) 

rim experiences a systematic drift, due to the probability "current" J; Jj — 
nmv{ij-i,ij,ij+i), arising from the deterministic evolution of each of the gaps ij at 
the rate of dij/dt = v{ij_i,ij,ij+i) = (1/2q;)(1/^^_i - 2/i^ + 1/i^+i), as given by 

— * 

Eq. (3) and absorbing a factor of 2 A in t. Then, using the notation im = {ii, ■ ■ ■ ,im), 
we can write the Fokker- Planck equation 

d ^ (9 r - 

j=2 i 



Q^J^ V{i',iui2)nm+l{i',im,t)di' 

v{im-U^rm^ )^m+l(^m)^ ,t)di 



d 



Jo 

m—1 



di 

m- 

- ^ hm [^;(£,-,£;,£,+i)n^+i(£i,...,£,-,£;,£,+i,...,Ct)] • (13) 

Note that since the expression for d£/dt requires knowledge of the gaps lengths to the left 
and right of i, the drift for ii (the second term on the rhs) forces us to consider m + 1 
gaps, including i' (the neighbour of ii to its left), and similarly for im (the third term 
on the rhs) which requires the gap i" to its right. The last term denotes the creation 
of nm when the gap i'^ shrinks and disappears (due to coalescence) from the (m + 1)- 
gap configuration £i, £2, • • • , ij-, i'ji ^j+i, ■ ■ ■ , ^m- Thus the Fokker-Planck equation for Um 
requires Um+i, resulting in an infinite hierarchy of equations. 

4.1. Scaling 

A useful simplification that gets rid of the time variable t results from the scaling 
assumption 

nm{il,...,im,t)=N{t)pm{il,...,im,t)=N{t)i-"'hm{Zu---,Zm); Zj=ij/i. (14) 
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In [1] we have presented simulation results supporting this scaling form for m = 1, or 
the single-gap pdf. For m = 2 we can test the scaling hypothesis in a roundabout way, 
by focusing on the distribution function of the particles' speeds 

/•OO /"OO 

p{v,t)= / 5[v~\^{l'-'^-t^n\)p2{l'A\t)M'df . (15) 
Jo Jo 

A simple change of variables now shows that (14) implies the scaling of the speed pdf: 
p{v,t) = i°'p{vi°'). In Fig. 2 we confirm that this scaling indeed holds, for the case of 
a = 1, giving us some confidence in the postulated scaling form (14) for general m. 



0.025 




Figure 2: Speed distribution function, shown for several times, for the case of a = 1. 
Inset: Collapse of the pdfs is shown for three different times after the scaling regime set 
in (by time T the initial number of particles has been reduced by a factor of 37), using 
the relation i ~ t^/^, valid for a = 1. 



Using (14) and the fact that N{t) ~ 1/^, the Fokker-Planck equation reads. 



d 



a{m + l)hmizi, . . . , Zm) + zj—hmizi, . . . 



(16) 



m—l r% 

y— 



i=2 
d_ 
dzi 



1/2 



dzj l\zf_i Zj 



1 1/2 
-r + 



hmi^Zi , . . . , Zjji^ 



^ (z' z z z )dz' 

a cT j '^"1+1 , Zl, Z2 ■ ■ ■ , Zm) CLZ 



2;'" z\ 



-2 / 
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dZr, 



JO 
m—1 

+ lim 



1/2 1 1/2 

— -\ — I ilm+l\Zi, 



/l/2_^ 



-la. 
^3 



1/2 



In fact, to arrive at this equation it is necessary to assume that £i°' is a constant (that 
can be absorbed in the time units, so we take it to be 1). This condition is in agreement 
with the coarsening law found from dimensional analysis, Eq. (8). 



4-2. Truncation scheme and Mean Field solution 

Eq. (16) is valid for m > 2. This infinite hierarchy of equations can be truncated at any 
level m by the usual Kirkwood approximation scheme: 

hm{zi, ■ ■ ■ , Zm)hm{z2, ■ ■ ■ , -^m+l) 



i^m-l[Z2, . . . , Zm) 

For example, at the lowest level of m = 2 we can truncate the hierarchy (16) using 

h2{zi, Z2)h2{z2, Zs) 



hsizi, Z2, Zs, 



h{z2) 



(17) 



(18) 



This yields 

d d 
3ah2{zi, Z2) + azi—h2{zi, Z2) + cxz2^h2{zi, Z2) 



^^{h2{z„Z2) 



d 



dzo 



h2{Zl,Z2) 



1/2 
1/2 



r dz' , . /1/2 1 

/ —h2{zx,A + ^ 

Jo Z \ Z2 Zi 



KZ2) 



%h2{z2,z')+(^-^ 



+ lim 



h2{z,, Z')h2{z2, z') (1/2 



1 



fl9) 

h{z') \z'^ z"^ ' z^ ) ' ^ ' 

which is a closed equation for /12, since h{z) = hi{z) = h2{z, z') dz'. This equation is 
hard to solve, or even to work out its limiting behavior. We simplify further by making 
the mean-field assumption that adjacent gaps are not correlated, so that 

h2{zi, Z2) = h{zi)h{z2) . (20) 

Inserting this in (19) and integrating over Z2, we obtain 

1 \ 1 



2ah{z) + azh'{z) 



h{z) 



^z" / z"- 

where the prime denotes differentiation with respect to z. Thus, 



(21) 



(22) 



h{z) az + Z-" - {z-") 

A naive inspection of Eq. (22) in the hmit of large z leads to the erroneous 
conclusion that h{z) ~ 1/z^. The problem is with {z^") (which must be determined 
self-consistently): depending on its value, the denominator on the rhs might vanish. 
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yielding a singularity that must be taken into account. Proceeding more carefully, we 
multiply Eq. (21) by and integrate over z, from to oo, to yield 

2a {z^) + az''+^h{z)\f^Q - a{k + 1) {z'') 

= [z\{z-'') - z-'')h{z)]Zo - k (^-"> {z''-^) + k (^'=-'-"> . 

Putting A; = we obtain the small- 2; behavior of h{z): 

{az°' for a > 0, 

') for « < 0. 

For k 00, and assuming that h{z) has finite support so that the various moments 
exist and limfc^oo {z^^^) / i^z^) = 1, we get, for the unknown (— aj-moment, 

(^-") = l + a;. (24) 

Armed with this information we now examine the denominator on the rhs of (22): For 
a > 0, az + z""' achieves a minimum of 1 + a, at z = Zc = 1, so that the denominator 
is always positive except at z = Zc, where it vanishes. Meanwhile, the numerator, 
—2a + Q;2;~°~^ is positive at z < Zm = (1/2)^/'^^+"', and negative for z > Zm, (note 
Zm < Zc). Therefore, h{z) rises as az"" (for z near 0), reaching a maximum at z = z^, 
then drops back to zero as z approaches Zc- The finite support assumed a priori is 
caused by the singularity at Zc and is now seen to be justified. For a < (but a > — 1), 
az + z'"' achieves a maximum of 1 + «, at 2; = = 1, so that the denominator is always 
negative except at z — Zc, where it vanishes. The numerator is negative for z < Zm and 
positive for z > z^, so h{z) starts at |q;|/(1 — \a\) at 2; = 0, rising to a maximum at 
z — Zm and drops back to zero towards z = Zc = 1. 

For the case of a < one can derive from (22) the correction to h{0): 

h(.) . MO) (1 + ^) = ^ (1 + ^.l"l) , „ < , (25) 

where we have used (23) and (24) for the last equality. Numerically, it is impossible to 
separate between the finite h{0) and the non-analytic correction, as both result in an 
infinite slope as 2; — 0. Indeed, the exponent j3 measured for a = —1/2 in Fig. 1 was 
obtained assuming h{0) = 0, yielding only an effective value. 

Finally, it is possible to obtain explicit solutions for several values of a. For example, 
for the special cases of a = 1,2, and —1/2 we get 

h{z) = —^-.e-'/^^-'\ 0<z<l; a = l, (26) 

[l-zy 

^(^)- (l-.)J(i + 2.)W9 ^"^^^^"^^ O^-^l; " = 2> (27) 

h{z) = ^ g-2Vi/{i-^) ^ 0<z<l; a = -1/2 . (28) 

(1 - V^) 

In all these cases (and others) it is easy to confirm the predicted small z behaviour, 
the location of the maximum at Zm = (l/2)"'^/*^^'*'"'', the finite support to z < Zc = 1, 
and that indeed (-2^") = 1 + a. In Figs. 3 and 4 we compare these mean-field results 
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to numerical simulations. Fig. 4 additionally demonstrates that the scaling regime is 
achieved, numerically, also for a < 0, though this takes a lot longer than for a > 0. 
Reaching scaling from an initial random configuration proved impractical and we had 
to start the system with all the gaps nearly equal but for tiny perturbations. 

For a 3> 1 the forces on particles across narrow gaps is immense and unaffected 
by neighbouring particles, therefore one expects that the mean-field prediction for small 
gaps improves as a increases (notice indeed the increasing agreement in /3 as a — * oo, 
in Fig. 1). For a < 0, on the other hand, it is the force across unusually large gaps 
that dominates and is unaffected by neighbouring particles, so there we expect that 
mean-field predicts the tail of the distribution increasingly better as a approaches —1. 



1.8 




h{z) 




(27) (solid curves) to actual simulations 
2 (b). In both cases z has been rescaled 



Figure 3: Comparison of h{z) from Eqs. (26), 
(symbols), for the cases of a = 1 (a) and a = 
so that (z) = 1. 



5. The Role of Correlations 



Consider now correlations between adjacent gaps i and £'. We expect that the two 
gaps are strongly anti-correlated: as the particle separating the two gaps moves one 
gap grows on expense of the other. This anti-correlation is clearly visible in Fig. 5, 
where we compare the (unconstrained) gap distribution to the gap distribution of i 
when i' < 0.5 {£) and when i' > 1.5 {£), for the case of a = 1. 

Despite the intractable equation for h2 it is possible to make a precise statement 
on the effect of correlations. We start with the exact equation 

-n{i,t) = --J^ v{e,i,nns{i',i,f,t)de'di", (29) 

and use the scaling hypothesis (14) to yield 

2aM.) + «4M=)4f fM.=.)-|^$^. (30) 
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Figure 4: Gap distribution for a = —1/2. Simulations are shown for times T (triangles), 
1.5T (squares), and 2.5T (circles) after the scaling regime had been reached: by time 
T the particles had been decimated 667-fold. The solid curve represents the mean-field 
result from Eq. (28). For the sake of comparison, z has been rescaled to yield {z) = 1. 



Integrating over z we get 



z" 



(31) 



2=0 



where we have used the notation 

1\ = rdz^h^iz^ 
^T/z Jo ^1 Hz) ' 
for the conditional average of " subject to the constraint that the interval ii = Zil is 
adjacent to an interval of length £ = zl. It is also understood that the limit oi z oo 
needs changing to z ^ Zc if there is finite support. 

Suppose that a > 0. It is easy to realise that the first term does not contribute for 
2; = and the second term does not contribute for ^ = 00. Therefore, we get 



a = lim h{z) 



+ lim 



h{z) 



-0 z° 



(32) 



If one neglects correlations the term (^Zi does not depend on z, h{z) vanishes for 
large z, and we found h{z) ~ az"' for small z, which agrees with Eq. (23). 
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If we take correlations into account, we know from simulations that h{z) ~ 2;^, 
with P > a (Fig. 1), and the second term vanishes. Therefore, (^r")^ ^lust diverge as 
l/h{z) as z ^ 00. Such behavior is believable, in view of Fig. 5, where we see that 
the conditional gap distribution when the adjacent gap is unusually large changes the 
f5 exponent for the small- 2; behavior. The change in j3 could result in divergence of the 
(—a) -moment in the required way. 



1.2 
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Figure 5: Anti-correlation between adjacent gaps ^ and I' (for a = 1). The unconstrained 
gap distribution of ^ (solid line) is compared to the gap distribution when the adjacent 
gap is unusually large, I' > 1.5 {i), and unusually narrow, i' < 0.5 {£). 



The following simple example serves to illustrate how such behavior may arise from 
anti-correlations. Suppose that two adjacent gaps, x and y, are perfectly anti-correlated, 
such that their sum is fixed: 

h2ix,y) x{2-x)6{x + y-2), 0<x,y<2, 

the (5-function ensuring that x -\-y = 2. Integrating over x we get h{y) ~ y{2 — y) (note 
that the example is symmetric; a similar h{x) is obtained integrating over y). Therefore, 



1 h{x,y) 
X Hy) 



dx 



lx{2 



xyi2-y) 



5{x + y — 2) dx 



2-2/ 
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which diverges in the same way as l/h{y), as y approaches the upper support hmit of 

yc = 2. 

Let us now anaUse Eq. (31) for negative a. The terms evaluated at z = cx) (or 
z — Zc) vanish and the term h{z)z^"'^ vanishes at 2; = as well. Therefore, we have 

/i(0)lim/zfl\ =|a| = /iMp(0)(4"l\ , (33) 

2— >0 \ / z \ / MF 

where we have used Eq. (23) for the last equality. Thus, for a < 0, correlations revise 
the value of /i(0). 



5.1. Input 

For homogeneous random input at rate R, we must add to the rhs of Eq. (29) the terms 



R(^-in{e,t)+2 n{i',t)di' 



The first term accounts for losses of a gap of length i, which are proportional to R£ and 
to the probability for having the gap in the first place, while the second term accounts 
for the creation of gaps when a particle is deposited inside a gap i' > i, exactly at 
distance i from either edge. Because of the randomizing effect of the input we expect 
that correlations play a lesser role and the mean-field approximation does somewhat 
better in this case. 

We specialize to the steady-state, setting the time derivative to zero, and use the 
mean-field ansatz to obtain 



a dz 



^)'^)h{z) -zh{z) + 2l\{z')dz', (34) 

where now z = i/is, {is ~ (A/i?)i/(2+a))_ r^Yie small-^ behaviour is very similar to 
that without input; h{z) - a (z) z'' for a > 0, and h{z) ~ h{0){l + ^1"!/ (-z'"')), 
h{0) = \a\ (z) I (^l"l>, for a < 0, while for large z mean-field predicts h{z) ~ exp(— ciz^) 
for q; > 0, and }i{z) ~ exp(— 02-2;^"'"') for a < (ci = a and C2 = |q;|/(2 — |q;|) are 
constants). These predictions seem to be borne out by simulations, although it is hard 
to make firm statements for the large-2; tail. The large-^ tail of the distribution for 
a = 1 provides yet another difference from diffusion-limited coalescence, where in the 
case of input p(z) ~ exp(— c^;^/^) [12]. 

In Fig. 6 we compare simulations and numerical integration of the mean-field 
equation for several values of a. The results are quite impressive, suggesting that input 
kills correlations very effectively. 



6. Summary and Discussion 

In summary, we have analysed a class of deterministic reactions where nearest particles 
attract as l/£" and coalesce upon encounter. The coarsening law is not sufficient to 
distinguish between different reaction dynamics so we have focused instead on the 
distribution function p(£, t) for the gap i between nearest particles. We have developed 
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Figure 6: Gap distributions for (a) a = 1, (b) a = 2, and (c) a = —1/2 in the case of 
input. Simulations (symbols) are compared to numerical integration of the mean-field 
equation (34). In all cases z has been rescaled to yield (z) = 1. 

an exact Fokker-Planck hierarchy of equations for the probabilities pm for finding m 
adjacent gaps of given sizes, capitalizing on the fact that the pm provide a complete 
description of a one-dimensional particles system. Truncating the hierarchy at the lowest 
level, we have studied the gap distribution at the mean-field approximation, neglecting 
all correlations between adjacent gaps. While this provides a qualitative description of 
the gap distribution, important details, such as the (3 exponent describing the small 
gap behavior p{i, t) ~ i^, could not be captured in this way. The role of correlations, 
however, was fully elucidated and we have found exact constraints that the distribution 
functions ought to obey due to their effect. The mean- field approximation worked more 
satisfactorily for the case of random input, where the role of correlations is diminished, 
yielding the same (3 exponent observed in simulations. 

Perhaps the largest outstanding problem is how to proceed, analytically, beyond 
the lowest mean-field truncation level, to obtain a more satisfying description of the gap 
distribution function and of the exponent /3. Especially intriguing is the numerical result 
that P ~ 1/2 for a = 0: Could this be derived exactly, maybe by a totally different 
approach? Cracking that question could suggest a more direct way of estimating /3 for 
all a values and perhaps for other deterministic reaction schemes. 
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